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Chromomagnetic Dipole-Operator Corrections in B — y X s j at 0{(3ool 2 ) 
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We calculate the fermionic corrections to the photon-energy spectrum of B — >• X s -y which arise from 
the self-interference of the chromomagnetic dipole operator Q& at next-to-next-to-leading order by 
applying naive non-abelianization. The resulting 0(/3oa^) correction to the B — > X s ^y branching 
ratio amounts to a relative shift of +0.12% (+0.27%) for a photon-energy cut of 1.6 GeV (1.0 GeV). 
We also comment on the potential size of resummation and non-perturbative effects related to Qg. 

PACS numbers: 12.38. Bx, 13.20.He 



I. INTRODUCTION 

The inclusive radiative i?-meson decay B — > X s ^/ rep- 
resents the "standard candle" of quark-flavor physics. It 
tests the electroweak structure of the underlying theory 
and provides information on the couplings and masses of 
heavy virtual particles appearing as intermediate states 
in and beyond the Standard Model (SM). See [ij f° r a 
concise overview. 

The present experimental world average for a photon- 
energy cut of E 7 > Eq with Eq = 1.6 GeV in the B- 



meson rest frame reads @ 



B(B -> X sl 



'exp 



(3.55 + 0.24 + 0.09) • 10" 4 . (1) 



The quoted value includes various measurements from 
CLEO, BaBar, and Belle @ and has a total error of be- 
low 8%, which consists of a combined statistical and sys- 
tematic error as well as a systematic uncertainty due to 
the shape function. 

In order to make full use of the available data, the SM 
calculation of B — > X s ^/ should be performed with similar 
or better precision. This goal can only be achieved with 
dedicated calculations of next-to-next-to-leading order 
(NNLO) QCD effects in renormalization-group improved 
perturbation theory. Considerable effort has gone into 
such computations. The necessary two- and three- loop 
matching was performed in Q and Q, while the mixing 
at three and four loops was calculated in and [7j . The 
two-loop matrix element including bremsstrahlungs cor- 
rections of the photonic dipole operator Q7 was found 
in 0, confirmed in and extended to include the full 
charm-quark mass dependence in jToj . The three- loop 
matrix elements of the current-current operators Q\,2 
were derived in within the large-/3o approximation. A 
calculation that goes beyond this approximation employs 
an interpolation in the charm-quark mass [l2l ]. Con- 
tributions involving a massive quark-loop insertion into 
the gluon propagator of the three-loop Q\,2 matrix ele- 
ments are also known [l3T | . Calculations of other missing 
NNLO pieces, such as the (Q7,Qs) interference were re- 
cently completed [l4[ . Further details on the status of the 



NNLO corrections to the branching ratio of B —y X s j can 
be found in (l5| . 

Combining the results listed above, it was possible to 
obtain the first estimate of the B — > X s j branching ratio 
at NNLO. For E = 1.6 GeV the result of the improved 
SM evaluation is given by [lU, [l6[ 1 



X s7 )^ >1 ' 6GeV = (3.15 + 0.23) 



10" 



(2) 



where the individual uncertainties from non-perturbative 
corrections (5%), parametric dependences (3%), higher- 
order perturbative effects (3%), and the interpolation in 
the charm-quark mass (3%) have been added in quadra- 
ture to obtain the total error. More details on the phe- 
nomenological NNLO analysis including the list of input 
parameters can be found in [12] ]. A systematic study of 
hadronic effects that cannot be described using a local op- 
erator product expansion has been recently carried out in 
[l7l ] (see also [HI). This analysis puts the naive estimate 
of the size of non-local power corrections in fl2l . [l6j on 
firm theoretical grounds, and at the same time indicates 
that a further reduction of the theoretical uncertainty 
plaguing ^ below 5% would require a theoretical break- 
through. 

Besides the branching ratio also the B —> X s j photon- 
energy spectrum is of theoretical interest and phe- 
nomenological relevance [l{|. While close to the phys- 
ical endpoint E y = jjib/2 the spectrum is dominated by 
the (Qi,2,Q7) and (Q7, Qi) contributions, the (Qs,Qs) 
interference is numerically the most important one for 
Ey < 1.1 GeV, because it involves a soft singularity 1/E 1 
related to photon bremsstrahlung. The theoretical de- 
scription of the (Qsi Qs) part of the spectrum has a sim- 
ple, but important feature, that is associated with the 
photon having a hadronic substructure, and manifests 
itself in the appearance of collincar singularities in the 



1 Several NNLO corrections (see 

EHHEl and p ar % 0) 

that 

were calculated after the publication of Il2i . [TB| are not included 
in the central value of {TJ, but remain within the quoted pertur- 
bative higher-order uncertainty of 3%. 



2 



perturbative result of the fixed-order calculation. The 
leading contribution of the (Qs,Qs) interference to the 
photon-energy spectrum in b — > Af art 7 has been known 
for some time [20(. This contribution is suppressed by 
a single power of a s with the respect to the leading 
(Q71Q7) interference, and therefore is part of the next- 
to-leading order (NLO) corrections to the spectrum. 

The (Qi,2,Qi,2), (<9i,2,Qr)) and (Qt,Q&) corrections 
to the photon-energy spectrum were calculated in the 
large- (3q approximation, i.e., including terms of order 
final through naive non-abelianization [2~j| . already in 
|22|. However, in that work neither the (Qi^Qs) nor 
the (Qs,Q 8 ) interference was considered. In this article, 
we close this gap partly by calculating the corrections 
to the photon-energy spectrum originating from the self- 
interference of the chromomagnetic dipole operator Q% 
at O(0o a s)- A calculation of the (Qi,2,Q8) interference, 



that completes the 0(a 2 s ) calculation of the spectrum in 
the large- /?o approximation has recently been completed 
and will soon be published (23j . 

This article is organized as follows. In Sec.|TT]we pro- 
vide the analytic results of our calculation, while Scc. lIIII 
contains a brief description of the calculation itself. The 
numerical impact of the considered NNLO corrections on 
the branching ratio of B — > X s ^ is studied in Sec. lIVI We 
conclude in Sec.lVl 



II. ANALYTIC RESULTS 

At the B-meson mass scale Hb = 0(m,b) the partonic 
b —> AP art 7 cut rate can be expressed in terms of the 
charmless semileptonic total decay width as 



J 



6a c 



v t * s v tb 



8 

£ 



(3) 



r 



where a om = ct om (0) = 1/137.036, 14; arc the rele- 
vant Cabibbo-Kobayashi-Maskawa matrix elements, and 
Cf s {nb) denote the effective Wilson coefficients defined 
as in fl2| . 

In the following, we will present analytic formulas for 
the 0(/3oa 2 ) corrections to K SS (E ). This function de- 
scribes the self-interference of the chromomagnetic dipole 
operator 



Qs 



9 



167T 



m b (ii) (s L a^T% R ) G a 



(4) 



where g is the strong-coupling constant, mb((J.) denotes 
the running MS mass of the bottom quark, qL,R are 



left- and right-chiral quark fields, G^„ is the gluonic field 
strength tensor, and T a are the color generators normal- 
ized to Tr (T a T b ) = T F 5 ab with T F = 1/2. 

Including QCD corrections up to NNLO, the coefficient 
Kss(Eq) can be written (in a notation following closely 
the one adopted in (j^ ) as follows: 



K 88 (Eo) = £ 



In agreement with |2C 
analytic expression 



47T 



|2J, we find for K$(E ) the 



(5) 



J 



4 

27 



In 



2 

mi 



5(2 + S) + 4 In 8] +4Li 2 <5 



2tt2 
3 



6(2 + 6) ln(5 + 81n<5 



26 3 



+ 3S 2 + 75 



(6) 



r 



where 5=1 — 2Eo/rrib and 8 = 1 — 6 = 2Eo/nib. As 

expected, the NLO function Kg$(Eo) is logarithmically 
divergent for both 8 — > (soft singularity) as well as 



m s — > (collincar singularity). Notice that terms sup- 
pressed by positive powers of the ratio m s /m,b have been 
neglected in ©. 



r 



(2) 

The NNLO function K y S8 '(E ) receives both fcrmionic 
and purely gluonic contributions. The former corrections 



arise from the Feynman diagrams shown in Figs. [TJ [21 
and [3l Since in the large-/3o approximation one con- 






FIG. 1: Four-particle cuts of the irreducible bottom-quark self-energy diagrams with a quark bubble contributing to the 
b — > sjqq (q = u,d,s) transition at 0(a 2 s ). Left-right reflected diagrams are not shown. The second and third diagrams give 
rise to collinear logarithms In (m%/rn%). 



siders exclusively massless fermion-loop insertions in the 
gluon propagators of the lower-order diagrams (2lj . it 
follows that in this approximation only the light-quark 
(q = u,d,s) loop diagrams shown in Fig.Q] need to be 
calculated. The graphs in Fig. [21 which also involve a 
light-quark bubble, belong to a new channel which opens 
up at NNLO, and therefore are not captured by the 
large-/3o approximation. Also notice that massless quark 
loops (q = u, d, s) in Fig. [3] involve a scaleless integral, 
which implies that they evaluate to zero in dimensional 
rcgularization. Consequently, we apply the naive non- 
abelianization prescription to the light-quark diagrams 

in Fig. [1] only and split Kg$(Eo) into a large-/?o and a 
remaining part 



^(2,/3 ) 
^■88 



(E ) + K, 



(2, rem) 
88 



(Eo) 



(7) 



and fa = 11 - 2/3 (N L + 2). As in the work [lj, we will 
set Nl = 3 in our numerical analysis. Notice that effects 
related to the absence of real charm-quark pair produc- 
tion in the partonic b — > Af art 7 decay and to non-zero 
values of the charm- and bottom-quark mass in quark 
bubbles on the gluon propagators are by definition con- 



tained in K 



(2, rem) 



(E ) and not in K^°\Eq). 



The coefficient kH' Nl \e ) 



L gg yj^u j introduced in ([8]) describes 
the contribution of the graphs in Fig.Q] involving a single 
massless quark flavor. It can be written as 



K, 



88 



(E ) 



dzF^-K^(E )L b 



(9) 



with 



A', 



(2,/3 ; 

88 



\P K^\E Q ). 



(8) 



where Lb = In f/^/m 2 ) 



2Ej/mb, and 



27 



in-:. 



36 - 41z + Viz 1 + 2z J 
6z 



lni 



1 + z 2 




H 


In 2 z 


z 


3 



60 - 65z + 16z 2 + 8z 3 , 604 - 702z + 126z 2 + 107z 3 
In z H 



6r 



36z 



(10) 



with z = 1 — z. The latter expression is the main analytic 
result of our paper. Similar to ^ also ([TU]) contains a 
collinear divergence, which we have regulated by keeping 
a non-vanishing strange-quark mass. Again terms sup- 
pressed by positive powers of m s /mb have been neglected 



in the function F> 



The function K 



,(2, rem) 



(E ) entering ([7]) encodes the 
0(a 2 ) contributions to the (Q$, Q$) interference that are 
beyond the large-/?o approximation. It takes the form 



J 



K 



(2, rem) 



(Eo) 



33 (2,N L ) 
TT -"-88 



(Eo)+ £ 



Kg^\E ) 



q—c,b 



K, 



(2,9,M) 
88 



(E a ) + Kg 9 \Eo) 



(11) 



The function K^ 9 \Eq) originates from diagrams with The calculation of these contributions is beyond the scope 
no quark loops, while K^ q,1 \Eo) corresponds to Fig. [2] 
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u,d,s 
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FIG. 2: Four-particle cuts of the irreducible bottom-quark 
self-energy diagrams with a light-quark bubble contributing 
to the b — } s^uu,dd,ss transition at 0{a 2 s ). Symmetric di- 
agrams are not shown. The shown diagrams give rise to 
collinear logarithms In (m^/m„ d s ) . In practice, these IR- 
sensitive terms are regulated by the light-meson masses M n 
and Mk ■ 



of the present article. As was already mentioned, real cc 



pair production is not included in b — > Xf^-f by defini- 
tion, while bb pair production is kinematically forbidden. 

Thus, K^' q (Eq) is non- vanishing for light quarks with 
q = u,d,s only. 

The function K 



(2,?,AQ 
88 



(Eq) originates from Fig. [3] and 
vanishes for massless quarks (q = u, d, s) in dimensional 
regularization due to the appearance of scalclcss inte- 
grals. Its analytic form for q = c,b can be obtained 
by multiplying the NLO coefficient (Eq) by a renor- 
malized one-loop vacuum-polarization function at zero- 
momentum transfer. Explicitly we find 



- ^^ T F K S8 



{E )L q . 



(12) 



where L c = In (/i^/mJV Notice finally that the effects of 
the diagrams in Fig.|3] can also been taken into account 
through gluon wave- function rcnormalization in the NLO 
graphs. 

It is also straightforward to derive an expression for the 



function A, 



(2, rem) 



(Eq) in the large-m c limit. The latter 



enters the calculation of the B —> X s -y branching ratio 
via an interpolation in the charm-quark mass [12j. In 
agreement with that paper, we obtain the expression 



J 



Kg^\E ) = ( -| + ^ - \l c ) Kg\E ) + Xg—\E ) . (13) 



Here the first term on the right-hand side is the lead- 
ing term in the large- m c expansion of Agg' rem ^(Ao). It 
consists of two parts, one arising from the normalization 
to the charmless semilcptonic rate and the other being 
proportional to L c , which is due to the MS matching 
corrections connecting the strong-coupling constants of 
the effective four- and five-flavor theories, as encoded 
in (|12l) . The m c -independent quantity xg' rem \E ) 
summarizes unknown 0(a 2 s ) contributions arising from 
the self- interference (Qs,Qs) in the theory with de- 
coupled charm quark (together with the corresponding 
bremsstrahlung) . 



III. CALCULATIONAL TECHNIQUE 

In order to calculate the O(f3oa 2 ) corrections to the 
partonic b — > A"P art 7 cut rate arising from the self- 
interference [Qs, Qs) we have employed the optical theo- 
rem. In particular, we have exploited the one-to-one cor- 
respondence between the interferences among diagrams 
contributing to the process b — ¥ sjqq and the physical 
cuts of three-loop bottom-quark self-energy diagrams. As 
can be seen by glancing at Fig.[TJ we are interested in di- 



agrams in which the chromomagnetic dipolc operator Qg 
appears on both sides of the cut. The contribution of a 
specific physical cut to the imaginary parts of the corre- 
sponding bottom-quark self-energy diagrams is thereby 
evaluated by means of the Cutkosky rules. See @ for 
more detailed discussions. 

We have evaluated the relevant four-particle cuts in 
two different ways. First, by a direct computation of 
the light-quark contributions using the set-up previously 
employed in the NNLO calculation of the (Q7, Q7.8) con- 
tributions [l, (3] , and, second, by performing the NLO 
calculation of the (Qs, Qs) contribution with a fictitious 
gluon mass which allows us to obtain the sought 0(a 2 s ) 
contributions from a dispersion integral over the gluon 
virtuality [2]|. 2 For a recent detailed review of this tech- 
nique we refer to (26[. In both cases, the reduction to 
master integrals via the Laporta algorithm [27| has been 
carried out keeping a non-vanishing strange-quark mass 



2 This method was also used in the calculation of the C(/3oc«s) 
corrections to the photon-energy spectrum of the (Qi,2iQi,2)i 
(Ql,2,Ql), and (Q7,Qs) terms [22l |. We verified the correctness 
of the (Qr, Qs) contribution given in the aforementioned article. 



FIG. 3: Three-particle cuts of the irreducible bottom-quark self-energy diagrams with a massive charm- and bottom-quark 
bubble contributing to the b — > sjg transition at 0(o? s ). Left-right reflected diagrams are not shown. The second and third 
diagrams give rise to collinear logarithms In (mj/mj). 



to regulate the residual collinear divergences. All the 
other infrared (IR) or ultraviolet divergences, appearing 
in intermediate stages of the calculation, have been regu- 
lated dimensionally in d — 4— 2 e dimensions. The master 
integrals have been evaluated analytically both by direct 
integration over the phase space and by employing the 



differential equation method [28j . Throughout the calcu- 
lation of the master integrals, we have neglected terms 
suppressed by positive powers of the ratio m s /mb- The 
agreement of the results obtained by the two methods 
serves as a powerful check of our calculation. 



J 



IV. NUMERICAL ANALYSIS 

In the following, we will investigate the numerical size 
of the 0{a 2 ) contributions to the {Qs,Qs) interference 
at the level of the branching ratio of B — > X s j. In order 
to simplify the comparison with the existing literature, 
we will adopt the conventions, the notations, and the nu- 
merical input parameters employed in Specifically, 
we will use (it, = 2.5 GcV, m& — 4.68 GeV, mb/m s = 50, 



B[B -> XJv) exp = 10.61%, C = 0.58, \V t * s V tb /V cb \ 2 = 

0.9676, a s (2.5 GeV) = 0.271, and Cg ff(0) (2.5 GcV) = 
—0.171. With this choice of input, one exactly repro- 
duces the central value of the SM prediction ([2]). 

We start by considering the impact of the large- /3q cor- 
rections. In this limit, we can write the correction to the 
B — > V s 7 branching ratio arising from the (Qs,Qs) in- 
terference at 0(a 2 ) as 



J 



AB{B -+ X sl ) 



= B(B -> X c lv\ 



6a c 



ttC 



v t * s v tb 



,ff(0) 



(Mb) I 



47T 



K, 



88 



(Eo), 



(14) 



r 



here C is the so-called semileptonic phase-space factor 
and K^'^°\Eq) has already been defined in ((5J). In the 
left panel of Fig.H we show AB(B -4 X s ^) E ~< >Eo nor- 



E^>E a 
SM 



as a 



malized to the central value of B(B 
function of the photon-energy cut Eq. We see from the 
solid line that the inclusion of the O((3 a 2 ) contributions 
leads to a relative change of the B — >• X s j branching ratio 
of +0.12% (+0.27%) for E a = 1.6 GcV (E = 1.0 GeV). 
We recall that for the two chosen values of E , the shifts 
due to the NLO corrections involving (Qs,Qs) amount 
to +0.24% and +0.66%, respectively. These numbers 
imply that after naive non-abelianization the term (|9]) 
constitutes a correction of almost 50% with respect to 
the 0(a s ) contributions. 

As we have already mentioned, in the (Qs,Qs) in- 
terference also the corrections which are not part of 



the large- /3o approximation (such as the four-particle 
cut diagrams in Fig. [5]) involve collinear logarithms as- 
sociated with photon fragmentation of b — > sg. Suf- 
ficiently far away from the endpoint of the photon- 
energy spectrum, the resulting Ill-sensitive terms can be 
subtracted and absorbed into non-perturbative photon- 
fragmentation functions [24|. 3 which obey perturbative 
evolution (Dokshitzer-Gribov-Lipatov-Altarisi-Parisi or 
DGLAP) equations with non-pcrturbativc initial distri- 
butions to be extracted from experiment. 4 While a com- 



3 In fact, in 1241 only purely perturbative corrections are included. 

4 In the endpoint region the non-perturbative physics associated 
with the (Qg,Qs) interference is encoded in a complicated sub- 
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E Q [GeV] E Q [GeV] 

FIG. 4: Left panel: Relative changes in B(B —t X s 'y)g^{ >Eo due to the O{j3oa 2 a ) contributions to (Qs, Qs). The solid (dashed) 
line shows the fixed-order (resummed) result as a function of the photon-energy cut Eo- Right panel: Comparison of perturbative 
and non-perturbative corrections related to the self-interference (Qs,Qs)- The solid line indicates the relative shift in B(B — ► 
^ s T)sm >£ ° due t° a resummation of collinear effects, while the dashed and dotted lines illustrate the impact of the non- 
perturbative components of the photon-fragmentation functions assuming two different models of vector-meson dominance. 
See text for further details. 



plete calculation of collinear effects at 0(a 2 s ) is beyond 
the scope of the present article, we find it illustrative 
to study the issue of IR-sensitive contributions arising 
in (J6j) , ([9]), and (fl"0|) . From this exercise we expect to 



get an idea about the potential size of both resummation 
and non-perturbative effects associated with the {Qs,, Qs) 
self-interference. 



The resummation of the collinear logarithms appearing 
in the K^' Nl \e ) corrections is achieved by convolut- 
ing the hard function C^ L (x), that describes the pro- 
cess b — > sqq for fixed energy x of the strange quark, 
with the universal strange-quark-to-photon fragmenta- 
tion function -D s _> 7 (x). Explicitly, we find that the re- 
sult of the resummation of the collinear logarithm in (fTU)) 
takes the form 5 



88 



with 



2ir 

^cm 



8 / 10 



1 dx 



C^(x)D s ^(^), (15) 



1 + x 



, (16) 



where x = 1 — x and [l/x] + denotes the usual plus distri- 
bution. A factorization formula similar to the one given 
in (fFo)) can also be derived for the complete 0(a 2 s ) cor- 
rection to Kss{Eq) in the collinear limit. 



5 In the absence of QCD, the expression for the strangc-quark- 
to-photon fragmentation function is given by D s _> 7 (x) = 
a em Qs/(2ir) (1 + x 2 )/x ln(fi 2 / fi 2 s ) with Q s = -1/3, fib ~ rn b , 
and fi s ~ m s « Aqcd- Substituting this result into (1151) . one 
recovers the terms in d 1 1) that are singular in the limit m s — > 0. 



The full photon- fragmentation functions -Di_>. 7 (x) with 
i = s,g are sums of perturbative and non-perturbative 
components. While the former are fully calculable in 
QCD, the latter have to be modeled. Following [29} . 
which the interested reader should consult for further 
details, we will employ a vector-meson dominance model 
and assume that quarks and gluons first fragment into 
vector mesons which then turn into photons. We be- 
gin our discussion by studying the impact of the anoma- 
lous parts of the photon- fragmentation functions, i.e., the 
components encoding the perturbative evolution as de- 
scribed by the inhomogencous DGLAP equations. Com- 
paring the resummed with the fixed-order 0(/3oa 2 ) re- 
sult, as indicated by the dashed and solid lines in the 
left panel of Fig. 01 respectively, we infer that the resum- 
mation of collinear logarithms decreases the obtained re- 
sults. Numerically, we find a relative change of +0.05% 
(+0.12%) for E = 1.6 GeV (E = 1.0 GeV), which im- 
plies that the resummation suppresses the considered cor- 
rection by more than a factor of 2. We also mention that 
for photon-energy cuts around 1.6 GeV the resummation 
of collinear logarithms appearing in the O^oa 2 .) correc- 
tion can be effectively described by choosing mi,/m s = 14 
in the analytic expression (fit)]) . 

We now turn our attention to the non-perturbative 
contributions related to the photon fragmentation from 
b — > sg. These corrections turn out to be potentially 
larger than the resummation effects. This is illustrated 



7 



by the right panel in Fig-H] which displays the relative 
change in B(B — > ^ s 7)sm >E ° arising from the sum of 
the 0(a s ) and O(/3 a%) corrections to (Qs, Q s ), including 
the anomalous parts of _D;_> 7 (a;) only (solid line) and em- 
ploying the full photon-fragmentation functions with two 
different non-pcrturbativc initial conditions (dashed and 
dotted lines). In each case, we have subtracted the fixed- 
order 0(a s ) corrections (|6]) from our results, since these 
effects are already part of the SM prediction ([2]). We 
see again that the choice mb/m s = 50, adopted through- 
out the recent literature on B — > X s j, tends to overesti- 
mate the effects of rcsumming collincar logarithms. Nu- 
merically, we find relative shifts of —0.07% and —0.23% 
for E = 1.6 GeV and E = 1.0 GeV, respectively. Af- 
ter incorporating on top of the anomalous also the non- 
perturbative components of Di^ f ~ j (x), we obtain instead 
corrections of -0.05% and -0.04% or -0.04% and 0.37%. 
The former (latter) numbers correspond to set I (II) of 
the full photon- fragmentation functions Dj_> 7 (x) deter- 
mined in [29j. We recall that while the initial condi- 
tions of the quark-to-photon fragmentation functions are 
well constrained by e + e~ data, the one of the gluon-to- 
photon fragmentation function is not. Compared to set 
I, the gluon-to-photon fragmentation function of set II is 
significantly larger, in particular, for small x. Since the 
function D g ^ y (x) enters the resummation of collinear 
logarithms in (Q S ,Q S ) at 0(a s ) via [13] 

K&\Eo) = — fdz \ [D s ^{z) + D g ^(z)} , (17) 

this results in larger shifts for set II than for set I. 

In conclusion, our study of non-perturbative effects re- 
lated to photon fragmentation seems to indicate, first, 
that the size of hadronic effects associated to the interfer- 
ence of (Qs,Qs) should not shift the central value of J2]) 
by more than +1% and, second, that setting mb/m s = 50 
in the terms In (m^/m 2 *) entering the fixed-order result 
allows one to capture most of the numerical effect. A 
recent much more detailed study [l?} finds slightly larger 
non-perturbative effects of [—0.3, +1.9]% related to the 
self-interference of the chromomagnetic dipole operator 
Qs- While a straightforward comparison of this result 
with ours is difficult, given the very different nature of 
the used framework, the fact that the two calculations 
result in numbers in the same ballpark gives us further 
confidence that hadronic contributions in B — > X s j re- 
lated to (Qs, Qs) indeed represent a minor effect. 

V. CONCLUSIONS 

In this work we have calculated the NNLO corrections 
to the b — > XP art 7 photon-energy spectrum in the large- 



st) approximation that arise from self-interference contri- 
bution of the chromomagnetic dipole operator Qs- The 
contributions from (Qs, Qs) ar e known to be numerically 
subleading at NLO for the photon-energy cut currently 
employed in the measurements of the B —> X s j branch- 
ing ratio. We find that this trend continues at NNLO 
and that the calculated O((3oa 2 ) corrections have only 
a marginal impact on the B X s j branching ratio, 
amounting to a relative shift of a few permille. However, 
corrections to the spectrum arising from the (Qs, Qs) in- 
terference are theoretically interesting in their own right, 
since they are known to be logarithmically divergent in 
the limit of vanishing photon energy, and because they 
contain collinear singularities that are associated with the 
intrinsic hadronic component of the photon. Concerning 
the latter issue, we have shown that non-perturbative 
effects in (Qs,Qs) due to photon fragmentation from 
b — > sg presumably constitute an effect of below a per- 
cent only. Our results can be readily incorporated in the 
SM calculation of the B — > X s ^ branching ratio. While 
a total non-perturbative uncertainty of about 5% will af- 
fect the SM prediction for the branching ratio for some 
time to come, it is still mandatory to update the avail- 
able NNLO estimate by including all the 0(a 2 s ) correc- 
tions which were calculated in the past four years, with 
the aim of reducing as much as possible the residual per- 
turbative uncertainty. The calculation presented here, 
constitutes a necessary ingredient to achieve this goal. 
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